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Motivated by the on-going rotating cryostat experiments in ISSP, Univ. of Tokyo, we ex- 
plore the textures and vortices in superfluid 3 He- A phase confined in narrow cylinders, whose 
qq ' radii are i?=50/Ltm and 115 /im. The calculations are based on the Ginzburg-Landau (GL) 

framework, which fully takes into account the orbital (/-vector) and spin (d-vector) degrees of 
freedom for chiral p-wave pairing superfluid. The GL free energy functional is solved numeri- 
q ' cally by using best known GL parameters appropriate for the actual experimental situations 

!— i ' at P=3.2MPa and i/=21.6mT . We identify the ground state /-vector configuration as radial 

^ ' disgyration (RD) texture with the polar core both at rest and low rotations and associated 

d- vector textures for both narrow cylinder systems under high magnetic fields. The RD which 
has a singularity at center, changes into Mermin-Ho texture above the critical rotation speed 
which is determined precisely, providing an experimental check for own proposal. 

■ 
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l o j . 

1. Introduction 

There has been much interest focused on superfluidity in various systems, ranging from 
superfluid 3 He liquid, 4 He super-solid, 1 - 1 neutral gases with bosonic ( 23 Na, 87 Rb, etc) and 
■ Fermionic atoms ( 6 Li), 2 ) to superconductivity with charged electrons, color superconductivity 

in dense quark matter. 3 ) Among them, superfluid 3 He 4 ~ 9 ) characterized by a triplet pairing 
occupied a special position because the rich internal degrees of freedom of a Cooper pair can 

> ■ 

be controlled by external nobs, such as the boundary condition of a confining wall, magnetic 
field and rotation. These features of He have been demonstrated experimentally and theo- 

a ~ 

retically, providing us a testing ground whose theoretical basis is well established to check 
various novel ideas and challenging proposals. Thorough studies in 3 He as a prototype multi- 
component order parameter (OP) superfluid are quite useful to understand other non-trivial 
pairing states, such as p-wave pairing superfluids via a Feshbach resonance, which is yet to 
be realized in ultracold atom gases. 10 ) Possible application of this study may be to a heavy 
Fermion superconductor UPt3, which consists of three superconducting phases in field H and 
temperature T plane. 11 ) This phase diagram can not obviously be explained in terms of a 
single OP, but must be multi-component OP. The high field and low T phase, so-called C 
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phase is analogous to the present A-phase. 

Superfluid 3 He-ABM(A) phase stabilized at high pressures and high temperatures (T) 
over BW(B) phase can be described by the ^-vector for the orbital p-wave symmetry and the 
(i-vector for the spin symmetry. These two vectors fully characterize the spatial structure, i.e. 
texture of the underlying OP symmetry. 4 ~ 9 ) 

It is known that in the absence of field, the Mermin-Ho (MH) texture 12 ' 13 ) is stable 
under confined geometries at rest where spontaneous mass current flows along the boundary 
wall. This MH texture is so generic and stable topologically, thus characteristic to multi- 
component OP superfluids. It exists even in the spinor BEC. 14 ~ 16 ) However, as the confining 
system becomes smaller and comparable to their characteristic length scale ~ ^ (~ 10//m), the 
texture formation becomes difficult due to the kinetic energy penalty, leading to destabilization 
of MH. Similarly MH becomes unstable under magnetic fields whose order is Hd (~2mT). 17 ) 
It is not known theoretically and experimentally what exactly condition needed for that and 
what kind of texture is stabilized in such a situation. 

Here we are going to solve this problem in connection with the ongoing experiments 
performed in ISSP, Univ. of Tokyo. They use the two narrow cylinders with the half radii 
R = 50^m and 115//m filled with 3 He-A phase. These two kinds of cylinders are rotated up 
to the maximum rotation speed f2=11.5rad/sec under a field applied along the rotation axis z 
and pressure _P=3.2MPa. They monitor the NMR spectrum to characterize textures created 
in a sample. So far the following facts are found: 18 ~ 22 ) 

(1) At rest the un-identified texture is seen for both samples with R = 50/xm and 115^m. 
The texture, which has a characteristic resonance spectrum, persists down to the lowest tem- 
perature T/T c =0.7, below which the B phase starts to appear, from the onset temperature 
T c = 2.3mK. Thus the ground state in those narrow cylinders is this un-identified texture. 

(2) Upon increasing 0, this texture is eventually changed into the MH texture which was 
identified before for R = 115/Ltm by the same NMR experiment. 20 ) This critical rotation speed 
S7 C ~ 0.5rad/sec, which is identified as a sudden intensity change of the main peak in the 
resonance spectrum. 

(3) With further increasing Q the so-called continuous unlocked vortex (CUV) are identified 
for R = 115/im sample. The successive transitions from the MH to one CUV, two CUV, etc in 
high rotation regions are explained basically by Takagi 23 ) who solves the same GL functional 
as in the present paper. However, the calculations are assumed to be the A-phase in the whole 
region. Thus the following calculation is consistent with those by Takagi. 23 ) Those successive 
transitions are absent in R = 50/xm sample because the estimated critical Q for CUV exceeds 
the maximum rotation speed 11.5rad/sec in the rotating cryostat in ISSP. 

(4) The un-identified texture for R = 115 fim sample stable at rest and under low rotations 
exhibits a hysteretic behavior about ±fl rotation sense, 19 ) meaning that this texture with 
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Cl > differs from that with Cl < 0, namely, this texture has polarity. These textures with ±fi 
can not continuously deform to each other. This feature is completely absent for the texture 
for R = 50^m sample because the NMR spectra for ±fi are identical. 22 ) Thus we have to 
distinguish two kinds of the texture stable at rest and under low rotations for i?=50^m and 
115/um. 

Before going into the detailed calculations, we introduce here the possible textures to be 
examined in the following: The texture consists of the orbital and spin parts, each of which 
is characterized by the /-vector and ci-vector respectively. Thus the total order parameter, or 
the texture is fully characterized by the combination of the /-texture and ci-texture. In Fig. 1 
we schematically show the possible /-textures and ci-textures where in MH the /-vector and in 
ax-type the d-vector change smoothly and flare out towards the wall. In contrast, the RD 24 ) 
(rd-type) has a singularity at the center. The /-vector in PA 25 ) (the cZ-vector in hb-type) shows 
a hyperbolic like spatial structure. In the following we examine the four textures, MH-ax, MH- 
hb, RD-ax and RD-hb, which turn out to be relevant and compete energetically each other in 
the present situations. 

The arrangement of the paper is as follows: After giving the Ginzburg-Landau (GL) free 
energy functional, we set up the GL parameters appropriate for the present experimental 
situation in § 2. Here the boundary condition which is essential for the present narrow cylinders 
is examined. We explain also our numerics to evaluate various possible textures. In § 3 we list 
up possible /-vector and d-vector textures both at rest and under rotation. In § 4 we identify 
the ground state texture by comparing the GL free energy for various radii R of the cylinder, 
rotation speeds CI and magnetic fields H. We critically examine the on-going experiments at 
ISSP by using the rotating cryostat in the light of the above calculations in § 5. The final 
section is devoted to conclusion and summary. 

2. Formulation 

2. 1 Order parameter and GL functional 

The OP of superfluid 3 He is given by a symmetric 2x2 matrix 4 9 ) 

A(p) = iA^a^o-yPi, (n,i = x,y,z), (1) 

where a is the Pauli matrix, p is the unit vector of the momentum on the Fermi surface. The 
summations over repeated indices are implied. Thus superfluid 3 He is characterized rank-2 
tensor components A^i inherent in the p-w&ve pairing (L = 1) with spin S = 1, where [i and 
i denote cartesian coordinates of the spin and orbital spaces, respectively. 

In order to understand the stable texture at rest and lower rotations, we examine the 
standard GL functional, which is well established thanks to the intensive theoretical and 
experimental studies over thirty years. 4-9 ) Namely, we start with the following GL form written 
in terms of the tensor A forming OP of p-wave pairing. The most general GL functional density 
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Fig. 1. (Color online) The possible textures of the /-vectors (upper) and d- vectors (lower) are shown 
schematically. The Mermine-Ho (MH), radial disgyration (RD), Pan- Am (PA) are for the /-vector 
textures (upper panel). The axial type (ax), radial disgyration (rd) and hyperbolic type (hb) are 
for the d- vector textures (lower panel). The points in MH, PA and rd-type denote positions of 
singularity. 



/bulk f° r the bulk condensation energy up to fourth order is written as 

/bulk = ~ a A^ ti A fl i + ftiA^A^AvjAvj + P2A fli A f j i iA U jA U j-\-P^A fli A ui A fJi jA l ,j 

+/?4 A^ j A V iA U j A^j +/?5 A^A U i A v j A^ j, ( 2 ) 

which is invariant under spin and real space rotations in addition to the gauge invariance 
U(l)xSO( 5 )(3)xSO^(3)- The coefficient a(T) of the second order invariant is T-dependent 
as usual, and the fourth order terms have five invariants with coefficients (3j in general. The 
gradient energy consisting of the three independent terms is given by 

/ grad = KiidZA^iPiArt) + K2(d*A* flj )(djA fli ) + K 3 (d* A^A.j). (3) 

where di = V, — i(^p-)(fi x r)j with the angular velocity f2 || z, which is parallel to the 
cylinder long axis and fi > means the counter clock-wise rotation. In addition, there are 
the dipole and magnetic field energies: 

/dipole = 9d (a^pAw + A*^ V A U ^ — -A^A^J^j , (4) 

/field = dmH^A^H^Ajyi. (5) 

In the A-phase the spin and orbital parts of the OP are factorized, i.e. A^ = d^Ai, where the 
spin part is denoted by and the orbital part by Ai. The d is a unit vector. We can describe 
the ABM and polar state under this notation. Now, Eqs.(2), (3), (4) and (5) are rewritten as 

/bulk = -aA*Ai + f3 13 A*A*A j A j + p 2 45A*A*A i A j , (6) 

/grad = K^d^idid^Aj) + K 2 (d*d fJ ,A*)(djd ll A i ) + K^d^A^d^Aj), (7) 
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dpdv (A^A U + A V A^ — —A U A U 



(8) 



2 



/field = 5 mA*^ d-H , (9) 



where /?i 3 = Pi+fh, and /3 2 45 = #2 + 04 + 05 ■ 

The coefficients a, 0j, Kj, g m , and ga are determined by Thuneberg 26 ) and Kita 27 ) as 
follows. The weak-coupling theory gives 

^mu.n^j^L), (10) 



3 V TJ U V T c/ 
_ORWC _ oWC _ R WC _ oWC _ _aWC _ 

7C(3)iV(0) 

-Mi x -p 2 -p 3 -p 4 --p 5 (ii) 

Kl = K 2 = K 3 = 7 ^ N ^ h \f^K, (12) 
J 240(7rfe B T c ) 2 ' V ; 

where N(0) and are the density of states per spin and the Fermi velocity, respectively. The 

coefficients ao and K are estimated by Eqs. (10) and (12) by using the values of N(0), T c and 

vf which are determined experimentally by Grey wall 28 ) within the weak-coupling theory. It 

is known that strong-coupling corrections for 0j are needed to stabilize A-phase, and we use 

the (3j values estimated by Sauls and Serene. 29 ) The value of ga is 30 ) 

2 



9d = m v 



/ 1.1339 x 0.45TfV , x 

(jftAT(O)ln , (13) 

where and 7 denote the permeability of vacuum and the gyromagnetic ratio, respectively, 
and Tp is the Fermi temperature defined by Tp = 3n/4iV(0)fcB with the density n. Finally, 
g m is given within the weak-coupling expression by 

7C(3)iV(0)( 7 ft) 2 

9m ~ 48[(l + F ^k B T c r { } 
with Fq the Landau parameter taken from Wheatley. 31 ) 

We set all the GL parameters to correspond to the above experimental pressure P = 
3.2MPa, which are summarized as 

3.81 x 10 50 J _1 m- 3 , (15) 

3 

f3 1 = -3.75 x 10 99 J" 3 m- 3 , (16) 

02 = 6.65 x 10 99 J" 3 m- 3 , (17) 

3 = 6.56 x 10 99 J" 3 m" 3 , (18) 

04 = 5.99 x 10 99 J~ 3 m~ 3 , (19) 
= -8.53 x 10 99 J~ 3 m- 3 , (20) 

K = 4.19 x 10 34 J _1 m _1 , (21) 



N(0) 

a =~ 
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g d = 5.61 x 10 44 J _1 m- 3 , (22) 

g rn = 1-35 x 10 44 J _1 m- 3 (mT)- 2 , (23) 

^ = 9.51 x 10 7 m- 2 s. (24) 



Corresponding to these parameters, the dipole field is estimated as Hd = \J gd/g m ~ 2mT, 
which is a characteristic magnetic field where the dipole and magnetic field energies become 
same order. 

The stable texture can be found by minimizing total free energy, 

F = J d 3 rf(r) = J d 3 r(/ bulk + / grad + / dipole + / field ). (25) 

We have identified stationary solutions by numerically solving the variational equations: 
5f{r)/5d^(r) = 0, 6f(r)/6Ai(r) = in cylindrical systems, assuming the uniformity towards 
the z direction. The height of the cylinders used in ISSP experiments are ~5mm, which allows 
us to adopt this assumption. Thus we obtain the stable d-textures and /-textures. Namely, we 
solve the coupled GL equations in two dimensions. The /-vector is defined as 



A*A k 

h = -ieijk | A , 2 , (26) 



where e^-fc is totally antisymmetric tensor and |A| 2 = A*A{ is the amplitude of OP. The mass 
current is given by 

ji = ^-ImiA^Aj + A*V 3 A t + A*V 3 A 3 ). (27) 

Here we expand the OP in a basis of spherical harmonics Yi m (I = 1, m = —1, 0, 1) for ease to 
express the initial configuration of the /-textures; 

A(p) = A + p + + A po + A-p-, (28) 

where p± = =F (Px ± Wy) /V?, Po = Pz and A± = =F (A x T iA y ) / V%, A = A z . 

It is noted that the characteristic length associated with the dipole energy is given by ^ = 



y/K/gd, which is estimated to be an order of 10/um. In contrast, the usual coherence length is 
£ = yj Kjoi = £o/\A — T/T c , where the coherence length at zero temperature £o ~ 0.01/zm. 
These two length scales differ by three order of magnitude which causes great difficulty to 
handle the problem involved both scales simultaneously. This is indeed our problem. The 
radial disgyration (RD) has a phase singularity at the center where the ^-vector vanishes 
around £-scale region while MH has no singularity and ^-vector is non-vanishing everywhere 
whose spatial variation is characterized by In order to compare two energies, we need 
to handle two scales simultaneously. Since at the boundary /-vector is constrained such that 
it is always perpendicular to the wall, thus RD and MH exhibit a similar /-vector texture, 
differing only around the center of a cylinder. We evaluate possible /-textures RD and MH in 
combination with ci-textures; axial type (ax) and hyperbolic type (hb). The radial disgyration 
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of d-texture (rd) is neglected since it has a singularity where superfluidity is broken. Namely, 
we mainly examine four kinds of texture, RD-ax, RD-hb, MH-ax and MH-hb in addition to 
the Pan- Am (PA) only for smaller systems in this paper. Those textures, we believe, exhausts 
all relevant stable textures. There is no other texture known in literature. 5 ) 

2.2 Numerics and boundary condition 

The numerical computations have been done by using polar coordinates (r, 9). The radial 
direction r is discretized into 1000 meshes for R=50/xm system and 2300 meshes for R=115/im 
while the azimuthal angle 6 is discretized in 180 points. Thus the total lattice points are 
1000x180 and 2300x180 in (r, 9) coordinate system for R=50/im and R=115/im respectively. 
The average lattice spacing is an order of 5£o f° r both cases, which is fine enough to accurately 
describe a singular core in RD. Note that we are considering high temperature region, £(T = 
0.95T C ) ~ 4.5£o ~ 45nm while our lattice spacing is 50nm. The advantages of using (r, 9) 
coordinate system over the rectangular (x, y) system are (1) we can reduce the total lattice 
points, keeping the numerical accuracy and (2) it is easy to take into account the boundary 
condition at the wall where the Z-vector orients along the perpendicular direction to the wall, 
namely the radial direction r. The disadvantage is that we can not describe the Pan- Am (PA) 
type configuration for the Z-vector. However, PA is topologically similar to RD because as 
seen from Fig. 1 two singularities at the wall in PA merge into a singularity at the center 
in RD. Thus we may ignore PA and concentrate on RD only except for smaller radius cases 
(R < 20/mi) where we confirm the above statement (see subsection 3.3 for detail). 

On each lattice points the 9 variational parameters are assigned, coming from the complex 
variable A^i = x, y, z) and the real three dimensional vector d^fx = x, y, z). These 9 parame- 
ters are determined iteratively and self-consistently by solving the coupled GL equations. One 
solution is needed ~7 days for R=50//m system and ~20 days for R=115/Ltm system by using 
OpenMP programming on a XEON 8 core machine. The applied field H = 21.6mT is taken 
for all calculation except for subsection 4.2. 

3. Stable Textures 

In order to help identifying the possible texture realized in narrow cylinders (i?=50/xm 
and 115/um) both at rest and under rotation, we first examine the detailed spatial structures 
for each texture, which consists of the /-vector and d-vector before discussing the relative 
stability among them under actual experimental setups. Here we study mainly four types of 
the textures: RD-hb, RD-ax, MH-hb and MH-ax, which are relevant for later discussions on 
experiments. As a supplement, we also examine PA-hb for smaller sizes. 
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3.1 Radial disgyration (RD) 

The RD is characterized by having a singularity at the center where A±(r = 0) = and 
Ao(r = 0) / 0. Thus the vortex core is filled by Aq component, namely it is a polar core 
vortex. 32 ^ The associated d-vector texture could be either hb-type or ax-type. The /-vector 
texture can be obtained by starting with an initial configuration: 

A + (r,0) = ^ ta nhC-)e- ie , 

Mr,9) = ^, (29) 

A_(r,0) = ^tanh(^) e *, 

where is the amplitude of OP in the bulk region far away from the center and boundary. 
The combination of the winding number in RD is (w+, wo,w-) = (—1,0, 1) where w+, wq, w- 
are the winding numbers for A + , Aq, A- respectively. The A± components vary over £ distance 
from the center while the Aq component stays constant throughout the system. Since at the 
center only the Aq component is non-vanishing, the polar state is realized there as mentioned 
before. It is seen from Eq. (26) that in this RD form l z = because of \A + \ = \A-\. 

In Fig. 2 we show the results of the stable RD-hb textures for R=50//m, where Figs. 2(a) 
and 2(b) at rest and Figs. 2(c) and 2(d) at CI = 3rad/s. It is seen that /-vectors flare out 
from the center and point perpendicular to the wall (Fig. 2(a)). The d- vectors point almost to 
the horizontal direction (the x-axis), and curve near the wall (Fig. 2(b)), because the dipole 
interaction tends to align the d- vector parallel to the /-vector direction. We call it "hyperbolic" 
(hb). The /-vectors do not point completely to the radial direction, but are twisted slightly 
due to the dipole interaction. Since the /-vectors are strongly constrained by the boundary 
condition, they should be perpendicular to the wall so as to suppress the perpendicular motion 
of the Cooper pair. In other words, one of the two point nodes which locate to the /-vector 
direction should direct to the wall, thus saving the condensation energy loss. Note that there 
is no z-components for both / and d in RD-hb at rest. It will be soon shown that this RD 
texture is most stable at rest and lower rotations in narrow cylinders. 

Once the rotation is turned on, the /-vectors and ci-vectors now acquire the z components. 
It is seen from Fig. 2(c) that the /-vectors tend to point the negative z direction, recognized as 
the color changes and also seen from Fig. 2(e). That region is confined along the y-axis where 
the dipole interaction, which tend to / || d, is not effective compared with other regions. In Fig. 
2(e) the cross sections of l z component are displayed. It is clear that l z component is mainly 
induced along the y axis. Simultaneously, as shown in Fig. 2(d) the d- vectors point to ±z 
directions whose boundary occurs along the vertical direction, that is, the y-axis. As shown 
in Fig. 2(f) the variation of the /-vectors produces the in-plane mass current. The je{x,y) 
distributions in the x and y plane are anisotropic. Along the radial direction the jg changes 
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Fig. 2. (Color online) Stable textures and current for RD-hb under R = 50/im and T/T c = 0.95. (a) 
/-vector textures at rest, (b) (/-vector textures at rest, (c) /-vector textures at fl = 3rad/s. (d) 
d- vector textures at ft = 3rad/s. The arrows indicate the x-y component and color codes denote its 
z component, (e) The z component of the /-vector and (f) the current jg along the radial direction 
r from the center. The red lines are at rest and blue lines are ft — 3rad/s. The solid (dashed) lines 
are along x (y) axis. 

the sign, namely, it is the same direction as the external rotation in the central region while 
opposite in the outer region. 

As for RD-ax, which is displayed in Fig. 3. The /-vector and d-vector textures at rest 
are shown in Figs. 3(a) and 3(b), respectively. These textures are seen to be cylindrically 
symmetric. The (i-vectors around the center point to the negative z direction (Fig. 3(b)). The 
area is characterized by magnetic coherence length £/j = w K/g m H 2 , which is estimated to 
be an order of 1 /im in H = 21.6mT. This length scale is larger than £, but smaller than 
Through the dipole interaction the z component of the /-vector is induced in that area (see 
Fig. 3(a) and red line in Fig. 3(e)). The mass current flows even at rest exclusively around 
the center (see the red line in Fig. 3(f)). Under rotation the z component of the /-vector in 
RD-ax acquires more negative component in the whole region (Fig. 3(c) and blue line in Fig. 
3(e)), so that the mass current increases further (blue line in Fig. 3(f)). The d-texture almost 
remains unchanged (Fig. 3(d)). 

3.2 Mermin-Ho (ME) 

The MH texture does not include singularities, being the ABM state throughout the 
system. The /-vector is directed to the z axis at the center, avoiding singularities. An initial 
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Fig. 3. (Color online) Textures and current for RD-ax under R = 50/im and T/T c = 0.95. (a) l-vector 
textures at rest, (b) <i-vector textures at rest, (c) Z-vector textures at 57 = 3rad/s. (d) d-vector 
textures at = 3rad/s. The arrows indicate the x-y component and color codes denote its z 
component, (e) The z component of the Z-vector and (f) the current jg along the radial direction 
r from the center. The red lines are at rest and blue lines are at fi = 3rad/s. 



configuration for the ^-vector texture is 

A + (r,6) = ^[l + cos (3(r)], 

A (r,e) = sin f3(r)e w , (30) 

A-(r,e) = ^[l-cos(3(r)}e 2ie , 

where (3{r) = irr/2R varies linearly from at the center to tt/2 at the wall. Thus the /-vector 
points towards the z direction at the center and the radial direction at the wall. The d-vector 
texture could be either hb-type or ax-type. 

The MH-hb is shown in Fig. 4 both at rest and under rotation. At rest the z component of 
the /-vectors is non-vanishing around the center in an anisotropic manner (as seen from Fig. 
4(a)). The anisotropy of the l z component in the x-y plane is understood in terms of the dipole 
interaction, which is shown in Fig. 4(c). The correlation between the Z-texture and d-texture 
produces the anisotropic arrangements in each textural configuration. The elongated l z com- 
ponent along the y direction is already explained in the previous subsection. The length scales 
of the l z variation towards the x and y directions are characterized by ^ and R respectively. 
The /-vector and the associated <i-vector variations along the x axis at y = around the 
center are shown schematically in Fig. 5. It is seen that since the Z-vectors (red arrows) flare 
out from the center, the z component of the d-vector (blue arrows) changes the sign at x = 0. 
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Fig. 4. (Color online) Stable textures and current for MH-hb under R = 50/im and T/T c — 0.95. (a) 
/-vector and (b) d-vector textures at rest. The arrows indicate the x-y component and color codes 
denote its z component, (c) The z component of the /-vector and (d) the current jg along the 
radial direction r from the center. The red lines are at rest and blue lines are £1 = 3rad/s. The 
solid (dashed) lines are along x (y) axis. 



-► X 



Fig. 5. (Color online) Arrangements of the /-vectors (red arrows) and d-vectors (blue arrows) for 
the MH-hb are shown schematically. The three pairs of arrows show schematically the variations 
along the x axis around the center. Note that the angle between the d-vectors and the x axis is 
exaggerated. 



This configuration embedded in the hyperbolic d-texture is most advantageous by increasing 
the parallel portion (/ \\ d) due to the dipole interaction. Since this structure is characterized 
by the winding number combination (w^,Wq,W-) = (0,1,2), it yields the spontaneous mass 
current at rest (red lines in Fig. 4(d)). Thus MH has polarity which breaks the symmetry for 
±17 under rotation, and jg(x,y = 0) is larger than jg(x = 0,y) around the center, and this 
relation is reversed in the outer region. Under rotation the overall configuration of the /-vector 
and d- vector textures are not much changed compared with those at rest. The small changes 
of the l z component and jg are shown in Figs. 4(c) and 4(d). 
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Fig. 6. (Color online) Textures and current for MH-ax under R = 50/im and T/T c = 0.95. (a) /-vector 
and (b) d- vector textures at rest. The arrows indicate the x-y component and color codes denote 
its z component, (c) The z component of the /-vector and (d) the current jg along the radial 
direction r from the center. The red lines are at rest and blue lines are f2 = 3rad/s. 



As for MH-ax shown in Fig. 6, the overall textures for /-vector and <i-vector have cylindrical 
symmetry. Note that the winding number combination (0, 1,2) is the same as before. The l- 
vectors and d- vectors around the center whose length scales are characterized by ^ and 
respectively. These two length scales also appear in jg as clearly seen from Fig. 6(d), where 
the sharp rise corresponds to and the maximum position to The l z component in Fig. 
6(c) shows gradual changes characterized by The external rotation hardly change these 
features. 



3.3 Pan- Am texture 

The /-texture in Pan-Am (PA) is characterized by the two singularities at the wall as 
shown schematically in Fig. 1. Here we show the stable PA texture for smaller size system 
R=20/xm in Fig. 7. In order to stabilize the PA we start with the initial configuration given 
by 

A + (r,9) = — ( snl ct — i cos a) , 

A (r,9)=i^, (31) 
j4_ (r, 9) = (sin a + i cos a) , 
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Fig. 7. (Color online) (a) /-vector and (b) (i-vector textures at rest for PA-hb under R = 20/im and 
T/T c — 0.95. The arrows indicate the x-y component. 



with 



tan a(r, 9) 



r 2 sin 29 



R 2 + r 2 cos 29 ' 

where R is the radius of the system. The angle tana diverges at (r, 9) = (R, ±7r/2). At those 
positions the OP amplitudes A± must be zero so that those are avoided. Around those points 
the phase windings are (w + , Wq, w_) = (—1,0, 1). Therefore those are analogous to the RD in 
the phase structure because the singularity at the center in RD is regarded to be split into 
two in PA. 

We show here the stabilized PA structure in Fig. 7 where the /-vector and d-vector struc- 
tures are displayed. From Fig. 7(a) the /-vectors tend to curve near the wall to satisfy the 
boundary condition, which causes the curving of the d- vectors shown in Fig. 7(b), otherwise 
they are straight pointing to the x direction. We notice that the energy of this PA is higher 
than the previous four kinds of the textures mentioned above. This is also true for other 
systems with different i?'s. In the following discussions we disregard the PA, only considering 
the previous four textures. 

4. Size-Dependence and Magnetic Field Effects 

4-1 Size- dependence and the critical rotation speed 

We show the free energy of four textures, which are RD-hb (blue), RD-ax (green), MH-hb 
(red) and MH-ax (light blue), in Fig. 8 for R = 50/um, R = 115/im and intermediate radius 
systems at rest. 

The most stable texture for the R = 50/iin system is the RD-hb. The reasons are un- 
derstood by the following. First we compare the free energy of the RD-hb with the MH-hb. 
The RD-hb is advantageous from the view point of the dipole energy because the /-vectors 
lie in the x-y plane together with the d-vectors. The winding numbers of the RD-hb are 
(w+j Wq, W-.) = (—1, 0, 1), so that the superflow velocity due to the phase gradient of the plus 
and the minus components is canceled. In contrast, since the winding numbers of the MH-hb 
are (w+ ,wq, W-) = (0, 1, 2), the MH-hb has the spontaneous current comprising of the super- 
flow velocity. Consequently the MH-hb costs the gradient energy at rest. On the other hand, 
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Fig. 8. (Color online) Size dependences of the free energies for various textures at rest. For R—bOfim 
RD-hb (blue) is stable while for i?=115/Ltm RD-ax (green) becomes stable. The other lines show 
MH-hb (red) and MH-ax (light blue). 



the RD-hb costs the condensation energy because the polar core with the radius £ exists at 
the center of the cylinder. Nevertheless, since the polar core is small enough compared with 
the system size, the RD-hb is more stable than the MH-hb. 

Next we compare the free energy of the RD-hb with the RD-ax. The gradient energy 
is favorable for the RD-hb since the d-vectors are almost uniform except near the wall. In 
addition, the magnetic field energy is gained because the cf-vectors lie in the x-y plane on 
application of the magnetic field toward z direction. On the other hand, the RD-ax is gained 
by the dipole energy because the (i-vectors are parallel to the /-vectors except with the radius 
£h around the center of the cylinder. Nevertheless, since the dipole-unlocked (the /-vectors do 
not parallel the (i-vectors) regions are small for the RD-hb, as well, the RD-hb is more stable 
than the RD-ax. For the same reason, the MH-hb is more stable than the MH-ax, and also 
the RD-hb is more stable than the MH-ax. Therefore the RD-hb is the most stable texture 
for the R = 50/^m system at rest. 

The most stable texture becomes the RD-ax for the R = 115/mi system. The reason is that 
the dipole-unlocked regions enlarge by increasing the system size for the RD-hb, whereas the 
regions of the RD-ax characterized by the magnetic coherence length £/, are almost constant. 
That is, the RD-ax is stabilized by the dipole energy for larger systems. Similarly the stability 
between the MH-ax and the MH-hb is interchanged. 

Under rotation, the MH textures are stable because they have the spontaneous current 
and they are tolerable under rotation compared with the RD textures which have very little 
spontaneous current. The external rotation drives the RD-hb into the MH-hb for the i?=50;um 
system. Similarly the most stable texture RD-ax for the R = 115/um system at rest changes 
into MH-hb by the rotation. Therefore /-texture and d-texture change at the critical rotation 
speed. 
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Fig. 9. (Color online) Critical rotation speed f2 c [rad/s] versus R[/im] between RD and MH. The 
dashed line is proportional to R~ 2 5 . RD is most stable than MH at rest and low rotations for 
various system sizes, including R—50fj,m and 115/im. 



The rotational speed changing the most stable Z-texture from the RD to the MH is defined 
as critical rotation speed O c . The Q c obtained from the numerical calculations for the systems 
with the various radii R is shown as red solid circles in Fig. 9. The RD is the most stable 
Z-texture from the R = 5^m to the R = 115//m systems at rest, and the MH-hb is the most 
stable texture under rotation. The VL C decreases with increasing the radius of the systems. 
The numerical results of the O c curve upward than the dashed line, which is proportional to 
i? -2,5 . Therefore it is concluded that the RD is the most stable texture up to the R —* oo 
systems at rest. Note that the results are obtained under the high magnetic field. It has been 
pointed out by Buchholtz and Fetter 17 ) that as the radius R of the systems is larger, the RD 
is more favorable under high magnetic fields. 

4-2 Magnetic field effect 

So far we have fixed the magnetic field at i?=21.6mT, which is selected by experiments 
at ISSP because the MNR signals are best resolved. From the theoretical point of view it is 
interesting to know in what weaker field MH is realized over RD. We show the result in Fig. 
10 where various textures are compared as a function of H for i?=50/um and T/T c =0.95. It 
is found that at H=2xn~£ MH-ax becomes stable over RD-hb beyond which RD-hb is always 
stable. Below the dipole magnetic field MH is expected to be stable over RD although it might 
be difficult to obtain good NMR signals in such a weak field. Thus for feasible magnetic field 
region around ff=21.6mT where sensible NMR signal is detectable, RD-hb is always stable 
in the narrow cylinders. 
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Fig. 10. (Color online) Comparison of various texture free energies as a function of the magnetic field 
H. The switching between RD-hb (blue) and MH-ax (light blue) occurs at H=2mT. The other 
lines show RD-ax (green) and MH-hb (red). 



5. Analysis of Experiments 

5.1 R=50fim 

In Fig. 11 we plot the critical rotation speed £l c as a function of temperature T/T c for 
-R=50/im case. At rest RD-hb is always stable over MH-hb. Upon increasing £1, MH-hb becomes 
stable at O c , which is an increasing function of T/T c because the vortex core energy ~ lni?/£ 
becomes lower as approaching T c where £ is longer. The reason why MH-hb is advantageous 
over RD-hb under rotation is that MH-hb has the spontaneous current, in contrast, RD-hb 
has no spontaneous current. These characteristics are seen in the inset of Fig. 11. 

NMR spectroscopy can be used to identify different topological objects of rotating su- 
perfluid 3 He-A phase, 33 ) namely it can be used to identify texture in a narrow cylinder. The 
dipole-locked (Z || d) and d _L H regions occupying most of the condensates in the cylinder 
yields the main peak of NMR spectra. In the high field limit its frequency is given by 4-6 ) 

u 2 t =u? L + ul (32) 

where col = jH is the Larmor frequency and oji is the longitudinal resonance frequency 
of the 3 He-A phase. The regions where the d-vectors are dipole-unlocked and deviate from 
the perpendicular direction to the field, generate the small satellite peaks. The resonance 
frequency of a satellite is expressed in terms of a relative frequency shift R 2 , defined by the 
equation 

u 2 t = u\ + Rluf. (33) 

Therefore the texture in a narrow cylinder is identified by finding R 2 . 

In experiment of the NMR spectroscopy for various temperatures a satellite peak with 
R 2 ~ 0.8 is observed at rest. 21 ) The value of R 2 corresponds approximately to that by the 
RD-hb. 34 ) We also notice that since RD-ax, which has a small area generating satellite peaks 



16/22 




0.65 0.7 



0.75 



0.8 0.85 

T/To 



0.9 



0.95 



Fig. 11. (Color online) Phase diagram of critical rotation speed O c (boundary line) and T/T c (R = 
50/im). At rest RD-hb is always stable. The inset shows an example of the energy crossing between 
RD-hb (blue) and MH-hb (red) at T/T c = 0.7, showing D, c = 1.8rad/s. 



around the center of the cylinder with a radius gives i? 2 ~ 1.0, the realized texture is 
definitely RD-hb, not RD-ax. Upon increasing f2, RD-hb changes into MH-hb. This texture 
has been identified previously by Takagi 23 ) who has calculated R\ ~ 0.3 for MH-hb. 

The precise determination of the temperature dependence of f2 c (T) is under way experi- 
mentally. Thus it is unable to check our prediction shown in Fig. 11 at this time. We emphasize 
here, however, that the NMR spectra at low £1 and high f2 are distinctively different, therefore 
it is clear that the phase transition between two textures RD-hb and MH-hb occurs. 

Experimentally 22 ^ there is no hysteresis for ±fi, indicating that the texture at rest has no 
polarity to the direction of rotation. This is in agreement with our identification of RD-hb for 
i?=50 / um, which has no polarity. 

5.2 R=115nm 

As shown in previous section, for the system size ii=115^m the ground state texture 
at rest is RD-ax which is stabler than MH-hb. The critical rotation speed fi c =1.3rad/s at 
T/T c =0.95 is calculated as shown in Fig. 12. The value of Q c roughly coincides with the ISSP 
experiment. 20 -* 

Experimentally there is a large hysteresis behavior about ±S~i centered at f2=0, 19 ) which 
is contrasted with R = 50/im case mentioned above. This is understood because RD-ax has 
polarity where Z-vectors and d-vectors point to one of the directions ±z-axis around the center 
which breaks ±fi symmetry. 

As for the NMR spectrum at rest, they observe no distinctive satellite feature for 
-R=115/xm. 21 ) These facts do not contradict our identification of RD-ax, which is i? 2 ~ 1.0. 
Thus we expect no satellite feature for RD-ax. 
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Fig. 12. (Color online) Free energy comparison of various textures for R = 115/im. RD-ax (green) 
changes to MH-hb (red) at £l c — 1.3rad/s. The other lines show RD-hb (blue) and MH-ax (light 
blue). 



6. Conclusion and Summary 

6.1 Conclusion 

Let us discuss the experimental facts (l)-(4) introduced in § 1. 

(1) We have identified the unknown texture as RD-hb for 50/im sample and RD-ax for 115/im 
sample at rest and lower rotations by finding the most stable texture in those conditions. 

(2) Upon increasing the rotation, we found that RD changes into MH where MH was identified 
before by NMR resonance shape unambiguously 19 ^ for 115/xm sample under rotations. This 
is confirmed by the present calculation. Namely, we precisely determined the critical rotation 
speed il c from RD to MH for two samples. 

(3) Under further higher rotation speeds multiple MH texture is expected to appear, that is, 
CUV. The estimated Q ~ 16rad/sec for R = 50/xm 18 -* is too high to attain in the present 
rotation cryostat in ISSP, whose maximum speed is 11.5rad/s. 

(4) Our identifications of RD-hb for R = 50/mi and RD-ax for R = 115/um are perfectly 
matched with the experimental facts that in the former (latter) there does not exist (do exist) 
the hysteresis under the ±fi rotations, namely, the former texture RD-hb has no polarity and 
can change continuously under the reversal of the rotation sense. In contrast, RD-ax has a 
definite polarity because the <i-vectors point to the negative z direction for the counter clock- 
wise rotation, which never continuously change into the positive z direction when the rotation 
sense is reversed, leading to a hysteresis phenomenon. 

Because of these facts, we conclude that the un-identified texture at rest and lower rota- 
tions confined in narrow cylinders is RD texture. Physically RD texture with the polar core 
at the center confined in a narrow cylinder under fields becomes energetically advantageous 
over MH, which is characterized by having the A-phase everywhere. In order to confirm our 
identification, we point out several experiments: 

(A) The most important prediction is the critical rotation f2 c (T) shown in Fig. 11 for R = 
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50/xm. As for R = 115^m, Q. c = 1.3rad/s at T/T c = 0.95 shown in Fig. 12. This can be 
detected by the change of NMR spectrum because RD and MH exhibit different spectral 
features. 33 ^ 

(B) It might be quite interesting to control the texture by tuning the magnetic field H as shown 
in Fig. 10. By lowering H, MH becomes stable at rest. This controllability by H should be 
utilized to identify textures, and furthermore be used for realization of exotic and unexplored 
physics associated with multi-component superfluidity, such as Majorana particle in parallel 
plates. 35 ' 36 ^ 

(C) The /-vectors in RD-hb and MH-hb, which do not point perfectly to the radial direction, 
exhibit a distortion or in-plane twisting as seen from Figs. 2(a) and 4(a). This ultimately leads 
to the mass current along the rotation axis; the z-direction. This non-trivial bending current 
(oc V x /) should be tested in a future experiment. 

6.2 Supplementary discussion 

Here we discuss the similarity and difference between the present superfluid 3 He-A and 
p-wave superfluids in atomic gases. 10 ) In order to gain some perspectives of the present study, 
it might be useful to discuss those issues. The similarity is obvious because both are described 
by the p-wave pairing. Thus the same GL functional forms, which is derived purely by general 
symmetry arguments can be used for both cases. Depending upon the situations, we can 
restrict the OP space to effectively reduce either orbital degrees of freedom or spin degrees of 
freedom for a Cooper pair. Thus we can regard the p-wave superfluids in atomic gases produced 
by a field sweep Feshbach resonance as a superfluid 3 He-A without the spin degrees of freedom 
because magnetic field polarizes the atom spin, thus the Fermion atoms are spinless. 

One of the main difference between them lies in the boundary condition. 3 He is usually 
confined by a rigid wall of a bucket while the atomic gases are trapped by a harmonic po- 
tential. The former strongly constraints the /-vector always perpendicular to a wall because 
the perpendicular motion of a Cooper pair to the wall is suppressed. In the harmonic trap 
the boundary condition acts more gradually and gently. The /-vectors tend to align tangen- 
tially at the boundary because the two point nodes along the /-vector direction are effectively 
excluded from the condensate volume region in this configuration, saving the condensation 
energy. Thus the resulting /-vector textures in two systems greatly differ. Under rotation the 
vortices emerged are also distinctive. 10 ) 
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